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Abstract 

The finite momentum meson spectral function (MSF) in the pseudoscalar channel 
is evaluated, adopting for the fermionic propagators HTL expressions. The different 
contributions to the meson spectral functions are clearly displayed. Our analysis may 
be of relevance for lattice studies of MSF based so far on the Maximum Entropy 
Method. As a further step the correlation function along the (imaginary-) temporal 
direction is evaluated. 
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1 Introduction 



In a previous paper [1] we studied the zero momentum Meson Spectral Func- 
tions (MSFs) and the correlations along the (imaginary) temporal direction 
in the QGP phase. 

Our starting point was the calculation done in Ref. [2], where the above quan- 
tities were evaluated in the Hard Thermal Loop (HTL) approximation. 
The main result of Ref. [2] was the unraveling of the behavior displayed by 
the MSF at small values of the energy u. Here indeed, owing to the presence 
of the plasmino mode in the HTL quark propagator, divergences appear in 
the density of states giving rise to peaks in the MSF, referred to as Van Hove 
singularities. 

The possible experimental relevance of such singularities was analyzed in 
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Ref. [3], where the back-to-back dilepton production rate was evaluated in 
terms of the zero momentum MSF in the vector channel. 
Actually in Ref. [2] HTL resummed propagators, in principle valid only for 
the soft modes, were adopted over the whole range of momenta. 
Hence in our previous paper [1] we compared these results with those obtained 
with a more realistic treatment of the hard fermionic modes. We found that 
the qualitative behavior of the MSF for small values of iv, arising from the 
presence of two branches (normal quark mode and plasmino) plus a contin- 
uum part in the HTL quark spectral function, was mildly affected. 

In the present paper we extend the analysis to the case of finite spatial mo- 
mentum. Since the problem turns out to be numerically more involved and 
since, as above stated, the qualitative structure of the MSF is substantially 
unchanged when corrections to the HTL scheme are accounted for, here we 
use HTL resummed fermionic propagators over the whole range of integration 
in momentum space. 

The analytical study of the spatial mesonic correlations was first addressed, 
in the non interacting case, in Ref. [4]. 

The study of these correlations is of particular relevance in finite temperature 
lattice QCD. The large number of lattice sites available along the spatial direc- 
tions allows one to study the large distance behavior of the mesonic correlators. 
This displays an exponential decay governed by a screening mass. Results ob- 
tained on the lattice can be found, for example, in Refs. [7,8,9,10]. Concerning 
analytical approaches, effects of the interaction on the meson screening masses 
were analyzed in Refs. [5,6] within a dimensional reduction framework. 
The knowledge of the screening masses leads to the understanding of the na- 
ture of the excitations characterizing the deconfined phase of QCD at the 
different temperatures: weakly interacting quarks and gluons or strongly cor- 
related states. 

In this paper we confine ourselves to the study of finite momentum mesonic 
correlators, leaving the spatial correlations (along the z-axis) for future work. 

Indeed the full information on the momentum and temperature behavior of 
the mesonic excitations in the QGP phase is encoded in their spectral func- 
tion. Unfortunately the latter is not measured directly on the lattice, but has 
to be reconstructed from the temporal correlator, which is known for a rather 
small set of points. This is the reason why most works are confined to the 
evaluation of the screening masses, even if the information which can be ex- 
tracted from the latter is much poorer with respect to the one stemming from 
the MSF. Nevertheless some results for the finite momentum MSF (for light 
mesons) start to be available [11,12]. 
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Presently the MSFs are extracted from the temporal correlators through a 
technique called Maximum Entropy Method (MEM). For a review of this ap- 
proach we refer the reader to [13]. 

We notice that, beside the study of meson properties at finite temperature, 
the problem of extracting spectral densities from lattice euclidean correlators 
is encountered in the evaluation of other quantities of interest for QGP phe- 
nomenology like the thermal dilepton production rate [14], the soft photon 
emission and the electrical conductivity [15], the hydrodynamical transport 
coefficients (shear and bulk viscosity and heat conductivity) [16,17,18]. The 
link between euclidean correlators, spectral densities and transport coefficients 
is provided, in the framework of linear response theory, by the Kubo relations 
[19,20,21]. 

Lattice results for finite momentum MSF are presently getting available also in 
the heavy meson sector. In particular the finite momentum charmonium spec- 
tral function (describing a cc meson moving in the heat bath frame) is found 
to displays a peaked structure in the pseudoscalar (rj c ) and vector (J/^) chan- 
nels even at temperatures above T c [22]. 

Since a J/$ moving in the heat bath frame sees more energetic gluons than 
at rest, the occurrence of a well defined peak in its spectral function seems 
to point to a robust binding of such a state since the collisions with "blue- 
shifted" gluons do not suffice to melt it. 

We hope that our results can provide a useful starting point or testing ground 
for MEM studies of MSF. 

Our paper is organized as follows. 

We shortly review in Sec. 2 the basic definitions of the mesonic operators, 

correlators and spectral functions addressed in this work. 

In Sec. 3 we briefly discuss the results of the HTL approximation for the quark 

propagator. 

Sec. 4 represents the central part of our paper. We focus our attention on the 
pseudoscalar channel. In this Section we present our analytical and numerical 
results both for the MSF (Sees. 4.1 and 4.2) and for the temporal correlator 
(Sec. 4.3). 

Finally in Sec. 5 we summarize and discuss our findings. 



2 Finite temperature meson spectral functions 



Following the notation adopted in Ref. [2] we consider the following current 
operator, carrying the quantum numbers of a meson: 

J M (-ir, x) = q(-ir, x)T M q(-ir, x) , (1) 
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where Tm = 1, 7 s , 7^, 7 M 7 5 for the scalar, pseudoscalar, vector and pseudo 
vector channels, respectively. We next define the fluctuation operator J M as 

J M (-ir, x) = J M (-ir, x) - (J M (-ir, x)) , (2) 

the average being taken on the grand canonical ensemble. 
The chief quantity we address is the thermal meson 2 point correlation func- 
tion: 



XM{-ir, x) = (J M {-iT, x)J ] M (0, 0)) 

= ( JmHt, x)jl(0, 0)> - (J M (-tr, x))(f M (0, 0)> 

^|j^e-^ XM K,P), (3) 

with r G [0,(3 = l/T] and uo n = 2nvrT (n = 0, ±1, ±2 . . . ). 

It is convenient to use the following spectral representation for the meson 

propagator in momentum space: 



XmK,p) = - / aw- a M {UJ,p) = -lmxM{uJ + ir],p). (4) 

J tU>„ — UJ 71 



-00 



Hence it is possible to express the thermal meson propagator in a mixed 
representation through the spectral function a M . Indeed starting from the 
definition: 



^ +00 

G M (-ir,p) = - ]T e- iuJnT X M(iu n ,p) 



n=—oo 

+00 



and performing the sum over the Matsubara frequencies with a standard con- 
tour integration in the complex uj plane [23,24], one obtains [2]: 



+00 i r / o /om +00 

G M (-ir,p)= / du a M (uJ,p) — = / du a M (u, p)K(u, r). 

J smh(c<;p/2) J 

(5) 

G(—ir,p) and (x(c<j,p) are the quantities of major interest for our study. We 
shall evaluate them in the rest of this paper. 

Also the z-axis correlator, which, in conformity with lattice calculations, is 
defined as follows 

G{z) = j dr J dx ± XM{-ir,x ± ,z) , (6) 



4 



can be expressed in terms of the finite momentum MSF a(u,p). 
In fact from Eqs. (3) and (4) it follows that: 



-oo 



J Z7T 



— oo 
+00 



^jP* z r°° duJ ^iP^lPA ( 7 ) 

27T J-00 uj K ' 



-00 



Clearly, due to rotational invariance, xm = Xm(\p\) and the choice of the axis 
is absolutely arbitrary: focusing on the z-axis one does not loose any general- 
ity 

It is justified to conjecture the spatial correlations to be exponentially sup- 
pressed at large values of z, namely: 



G{Z) ~ e ~ m scrZ _ (g) 

2^+00 

Hence from the evaluation of the meson screening mass m scr , governing the 
large distance behavior of the correlator, one can extract informations on the 
nature of the excitations characterizing the QGP phase. 
In the rest of this paper we shall show how the approach employed in Refs. [1,2] 
can be extended to the evaluation of the finite momentum meson spectral 
functions and temporal correlators. As already anticipated, due to the hard 
numerical work required to get Q(z), the evaluation of the z-axis correlations 
is left out for future work. 



3 HTL quark spectral function 



Fig. 1. The diagram contributing to the quark self-energy in the HTL approximation. 
The wavy line denotes a gluon. The external quark lines are soft. 

In this Section, following [24], we recall the spectral representation of the 
HTL resummed quark propagator which will be exploited in the evaluation of 
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the mesonic correlators and will result helpful in understanding the physical 
meaning of the different terms contributing to the MSF. 

We remind the reader that the HTL approximation for the quark propagator 
amounts to consider the self-energy diagram in Fig. 1. 

It can be shown that the major role in dressing the quark propagator is not 
played by the vacuum fluctuations (a quark emitting and reabsorbing a gluon), 
but by the interaction with the other particles of the thermal bath (antiquarks 
and gluons), carrying typical (hard) momenta K ~ T. 
The HTL quark propagator can be cast in the following form [24] : 

*S(u>,p) = *A + (uj,p f° "J"* 3 + *A-(u,p) ' f+ ?' P ( 9 ) 



being 

*A±(u,p) = -J 



m 9 

2p 



l^)ln^±2 
pi oj-p 



(10) 



with the quark thermal mass m q = g (T)T / 'y/E , g(T) being the gauge running 
coupling evaluated at the renormalization scale /i~T. 
Likewise the following representation 

s^p) = -tzlil 7 d j^hA _ i!±zve 7 d jJ^A , (11) 

2 J iuj n — uo 2 J iuj n — uo 

— oo —oo 

is easily proved to hold as well. In (11) p± and *A ± are related as follows: 
*A±(z,p) = - { dcj^^- =>p±(u,p) = -Im*A ± (" + ^P) , (12) 

J Z — UJ 71 



-oc 



with z off the real axis. 

Lumping together the above results one gets the following compact expression 
for the HTL quark spectral function 

P htl(co,p) = 7 ° ~ 2 7 ' P P+ (^P) + 7 ° + 2 7 'V (^P) , (13) 



leading to the spectral representation 



+ 0O 



*Qf ^ [a Phtl(^,p) n/l x 
S^u; n ,p) = - / dw— , (14) 

J lUJ n — UJ 

— oo 

for the thermal fermionic Matsubara propagator. 

The explicit expression of the HTL quark spectral function introduced in Eq. 
(11) reads [23]: 

2 U2 

P±(oo, k) = ^-^-[S(u - u±) + 5(cu + uj t )] + p±(u, k)9(k 2 - uj 2 ) , (15) 
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Fig. 2. Dispersion relations corresponding to the quasiparticle poles of the HTL 
fermion propagator in the time- like domain. We use the dimensionless units k/m q 
and uj/m q . 




0.05 - 



(0/m 



m„ 



Fig. 3. Dimensionless spectral function m q p± for space-like momenta at p 
as a function of uj/m q . The maximum of (3- arises from a second solution ui-(p) 
of the equation Re(*Al 1 ) = in the space-like region, but in this case, due to 
the imaginary part of the self-energy, it does not give rise to a well defined peak 
corresponding to a quasi-particle excitation. 

with 



P±(u,k) = 



±u> — k 



k{-u±k)+m* (±l-^lng) 



+ 



IT 2 ±w — k 
2 ' n q k 



2 ' 



(16) 

where the two dispersion relations cu±(k) are shown in Fig. 2, while the be- 
havior of f3±(u>,k) is plotted in Fig. 3. The thermal gap mass of the quark 
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Fig. 4. Diagram giving the pseudo-scalar meson 2-point function in the HTL approx- 
imation. The interaction vertices are 7 5 matrices, while the dressed lines are HTL 
resummed propagators. No vertex correction is required within this scheme. 

m q = gT / v6 shows up in Fig. 2 at k = 0. 

The HTL spectral functions in Eq. (15) consist of two pieces: a pole term, 
arising from the zeros of the real part of the denominators of *A± in Eq. 
(9), and a continuum term, corresponding to the Landau damping of a quark 
propagating in the thermal bath. 

In the time-like domain uj > k the quark spectral function is characterized by 
the quasiparticle poles uj = u±(k). The dispersion relation uj + (k) satisfies the 
equation *A^_ 1 (uj + (k), k) = and corresponds to the propagation of a quasi- 
particle with chirality and helicity eigenvalues of the same sign. The dispersion 
relation u-(k) satisfies instead the equation *A: 1 (w_(^), k) — and describes 
the propagation of an excitation, referred to as plasmino, with negative helicity 
over chirality ratio. Note that both these excitations are undamped at this level 
of approximation, since the logarithm in *A± doesn't develop any imaginary 
part in the time-like domain. 

On the other hand, in the space-like domain uj < k, the quark spectral function 
gets the continuum contribution /3±(uj,k) (Landau damping). Formally such 
a contribution is related to the imaginary part developed by the logarithm 
contained in Eq. (10). Physically it reflects two different processes: a quark 
absorbing a gluon from the thermal bath and being scattered by the latter; a 
quark annihilating with an antiquark from the thermal bath to give a gluon. 
In the following we investigate how this structure of the finite temperature 
quark spectral function affects the finite momentum mesonic correlations. 



4 Mesonic correlations in the pseudoscalar channel 

In the following we shall confine ourselves to the pseudoscalar channel. 
As shown in [25], the pseudoscalar vertex receives no HTL correction . Hence 
the HTL approximation for the meson 2-point function in the pseudoscalar 
channel, graphically displayed in Fig. 4, amounts to simply employ for the 
fermionic lines the resummed propagators given by Eq. (9). This leads to the 
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following expression: 

1 + °° r d 3 k 

X ps (icu h p) = 2N c - J2 J ^Tr[ 7 5 ^(^ n ,fc) 7 5 ^ n - ? ^,fc-p)] . (17) 

Setting q = k — p and making use of the spectral representation (14) of the 
quark propagator one gets: 



i +oo j3i, +°° +°° 

x ^ hP) = 2N c -^j—J^j^ 



lUJ n — UJ\ lUJ n — lUJl — U>2 



■X 



x Tr[7 5 pHTL(^i, fc)7 5 pHTL(^ 2 , <?)] . (18) 



The above allows us to extend the work done in Ref. [1] for the MSF and the 
temporal correlator to the case of finite spatial momentum. 



4-1 Finite momentum spectral function: analytical expressions 



Summing over the Matsubara frequencies in Eq. (18) with a standard contour 
integration, performing the usual analytical continuation iuj\ — > uo + ii] + (re- 
tarded boundary conditions) and taking the imaginary part of the result thus 
obtained one gets for the MSF the expression: 



+oo +OC 



a^(u,p) = -2N c J— -^(e^-1) J du, J dw 2 n(c;i)ra(c; 2 )x 

— oo — oo 

x 6(u - 0J\ - cu 2 ) ■ Tr[7 5 pHTL(^i, fc)7 5 PHTL(-^2, q)] • (19) 
Now, inserting Eq. (13) into Eq. (19) and since 



Tr 
Tr 



' 5 7°T7 • fe 5 7° =F 7 • Q 



' 5 7°T7 • fc 5 7° ± 7 • Q 



-(1-k-q), 
-(1 + k-q), 



(20a) 
(20b) 



one gets [2] 



o- ps (u,p)=2N c J 



al 3 k 



H-oo +oo 



1) J dui j dLU 2 n(uJi)n(LU2)5(iu — u)\ — cu 2 ) x 



x {(1 + fc • q)[p+(u 1 ,k)p-(-u 2 ,q) + p-(u 1 ,k)p+(-u 2 ,q)}+ 

+ (1 - fe • q)\p + (u u k)p+(-U2,q) +p-(u u k)p-(-U2,q)]} . (21) 
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The above expression can also be written as 



a ps (iu,p)=2N c J-^-(e f3uJ - 1) J du x J dw 2 n{uj 1 )n{uj 2 )5{^ -w 2 ) x 

— oo — oo 

x {(1 + fe • q)[p + (u 1 ,k)p + (u 2 ,q) + p-(u 1 ,k)p_(u 2 ,q)} + 

+ q)\p + (u 1 ,k)p-(u 2 ,q)+P-(wi,k)p+(u 2 ,q)]} ■ (22) 



where the identity p+(— u, k) = p-(ou, k) has been used. 

As discussed in Sec. 3 the HTL quark spectral function p±(u,k), whose ex- 
pression in given in Eq. (15), reflects the singularities of the quark propagator 
(10) in the complex o;-plane. These lie on the real tu-axis. At a given value of 
the spatial momentum k, in the time-like domain {oo 2 > k 2 ) discrete poles are 
associated to quasiparticle excitations; a cut for space- like momenta (u 2 < k 2 ) 
is related to the Landau damping. 

Inserting the explicit expressions for p±(cu,k) into Eq. (21) one gets then 
three different contributions (pole-pole, pole-cut and cut-cut) to the HTL 
pseudoscalar spectral MSF, namely: 



< T >,p) = a^(u,p) + a^(co,p) + a cc (u,p) . (23) 



Below we give the explicit expressions for these three terms. 
In the following 



(24) 



are the residues of the quasi-particle poles. 

Moreover we shall denote with q + (q + ) the normal quark (antiquark) mode, 
with q_ (q_) the plasmino (antiplasmino) mode and with M the excitation 
carrying the quantum numbers of a pseudoscalar meson. 
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The pole-pole contribution (for u > 0) is made up of twelve terms and reads: 



ax 



1 

k\ 



2(1 + fei • qfi)n(w+(A;i))n(w + (gi))Z + (A;i)Z + (gi) 1 - F 



fei 



+ 2(1 + fc 2 • qf 2 )n(w+(A; 2 ))[l - n{oj_{q 2 ))}Z+{k 2 )Z_{q 2 ) — 



/bo 



k 2 



u' + (k 2 ) - u;'_(q 2 )\ 



+ 2(1 + k ■ q 3 )[l - h{u_{k 3 ))]n{u + {q 3 ))Z_{k 3 )Z + {q 3 ) — h 



k 3 

k 2 

+ Y / (i + k4-qdH^4h))n(uJ-{Qd)Z-(h)Z_{q 4 ) — 



tjj'-{k 3 ) - w' + (q 3 )\ 



k4 



u'_(h) +wl(g 4 )| 



+ 2(1 + fcs • g 5 )n(o;_(A; 5 ))[l - n(^ + (g 5 ))]Z_(A; 5 )Z + (g 5 )— h > 



*>5 



ijj'-{k b ) - uJ+(q 5 )\ 



+ 2(1 + k ■ q 6 )[l - h(u; + (k 6 ))}h(u;4q 6 ))Z + (k 6 )Z_(q 6 )— H 



k@ 



uJ+(k 6 ) - u'_(qe)\ 



+ 2(1 - k ■ q 7 )h{u + {k 7 ))h{u_{q 7 ))Z + {k 7 )Z_{q 7 ) — ^ 



k 7 



uj' + {k 7 ) + (v'_(q 7 )\ 



+ 2(1 - k • q 8 )n(cu + (k 8 ))[l - h(cu + (q 8 ))]Z + (k 8 )Z + (q 8 ) — ' ' 



k& 



u' + (k 8 ) -u' + (q 8 )\ 

+ 2(1 " k ■ q 9 )[l - n(u4k 9 ))]n(u.(g s ))Z.(k 9 )Z.(g 9 ) % 

k 9 \uj_{k 9 ) - u_{q 9 )\ 

k 2 

+ 2(1 - fcio • q w )n{oj-{k w ))n{oj + {q w ))Z_{k w )Z + {q w ) - , n r~~T7 — T\ + 



kio 



\u'-{k w ) +u)' + (q m ) 



2(l-fcn • g 11 )n(fa;-(A ;il ))[l-n(a;-( gll ))]Z-(fc 11 )Z-( gll ) f 11 + 

k 2 

2(1 - fci2 • ai 2 )[l - w(a; + (fci 2 ))]w(u; + (gi 2 ))Z + (fci 2 )Z + (gi2) | , r 12 , r 

(25) 



In the above the q n = yk 2 +p 2 — 2pk n x and the k n are then the solutions 
of the following equations, corresponding to the physical processes explicitly 
specified below: 
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1) uj — ou+(ki 

2) uj — uj + (k 2 

3) uj + uj~(k 3 

4) uj — cj„(A; 4 

5) uj — uJ-(k 5 

6) uj + uj + (k 6 

7) uj — uJ+(k 7 

8) uj — uj + (k 8 

9) uj + uj-(k 9 

10) uj — uJ^(k 1Q 

11) u; — w_(/cn 

12) uj + cu + (A;i2 



- = 

+ uj-(q 2 ) = 
-uj+(q 3 ) = 

- cj_(g 4 ) = 
+ u + (q 5 ) = 
-w_(g 6 ) = 

- cj_(g 7 ) = 
+ uj+(q s ) = 
-U-(q g ) = 

- w+(gi ) = 
+ uj-(q u ) = 

- uJ+iqu) = 



<?+ + q + -> M 
g+ -> g_ + M 
q+^q. + M 
g_ + g_ -> M 
g_ -> g+ + M 
g_ -> g+ + M 
q + + q^—>M 
q + ^q + + M 
q- —> q- + M 
q- + q+^ M 
g_ -> g_ + M 
g+ -> g+ + M . 



The above identifications follow quite easily after inserting into Eq. (21) the 
quasiparticle contributions of the quark spectral functions given by the equa- 
tions 



pP° le (o>, k) = Z + {k)5{uj - uj + {k)) + Z_{k)5{uj + w-(k)) (27a) 



quark mode antiplasmino 

pole / 



p p _ ole (uj, k) = Z_(k)S(u - u-(k)) + Z + (k)5(uj + u>+(k)) . (27b) 



plasmino antiquark mode 



The role played by the statistical distributions in Eq. (25) appears then clear. 
In the annihilation processes the initial particles belong to occupied states 
(hence the product of two Fermi distributions). In the decay processes the 
initial particle belongs to an occupied state (hence the Fermi distribution), 
while the final one is produced in an empty state (hence the Pauli blocking 
factor). 
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The pole-cut contribution is found to consist of sixteen terms and reads: 



N 



(w,p) = ^(e^ - 1) J dkk 2 J dcos6{(l + k-q)x 



o -1 

.2 /, , , , /u\\2\ 



x Z + (k)(3 + (uo - u+(k), q)n(u + (k))n{u - uj + (k))6(q 2 -(to- u + (k)) 2 ) + 

+ Z_{k)f3+{uo + u-{k), q)(l - h{u-{k)))n{u + uo-{k))d(q 2 - {u + cu_(A;)) 2 ) + 
+ Z+(q)P + (w - u + (q),k)n(uj + (q))n(u - u + (q))6(k 2 - (u - u + (q)) 2 ) + 

+ Z-(q)P+(u + u-(q), fc)(l - n(cu_(g)))n(cu + cj_(g))#(A; 2 - (u + u>-(q)) 2 )+ 
+ Z_(A;)/3_(^ - u-(k), q)n(u-(k))n(u - cu4k))9(q 2 - (u - u_{k)) 2 )+ 

+ Z + {k)P-{u + w+(A;), q)(l - n(w+(*0))"(" + cj+(A;))% 2 - (u + w+(A;)) 2 ) + 
Z_(q)P_(u - A;)n(^(g))n(^ - u_(q))9(k 2 - (u - u>-(q)) 2 )+ 

+Z+ (g)/L(c<; + fc)(l - n(w + (g)))n(w + co + (g))fl(A; 2 - (lu + co + (g)) 2 ) 

+ (1 -fc • g)x 



+ 



x 



Z + (k)P_(uj - co + (k), q)n(uj + (k))n(uj - uj + (k))6(q 2 - (w - uj + (k)) 2 ) + 

+ Z_{k)P-{oj + tv-(k), q)(l - n{uj-{k)))n{uj + u^k))9(q 2 - (u + oj-(k)) 2 )+ 
+ Z_(q)P + (u - cu_(q), h)n(u-(q))n(u - uj^(q))6(k 2 -(to- UJ_(q)) 2 )+ 

+ Z + (q)p + (uj + u+(q), fc)(l - n(u+(q)))n(u + uj + (q))9(k 2 - (u + co+(<?)) 2 ) + 
+ Z_{k)P + {u - w_(ifc), q)n(u-(k))n(u - uu_(k))d(q 2 - (u - u_(k)) 2 ) + 

+ Z + (k)P + (u + u + (k), q)(l - n{uj+{k)))n{uj + uj + (k))6(q 2 - (w + u+(k)) 2 ) + 
Z + (q)(3_(u - u+(q), k)n(uj + (q))n(uj - uj + (q))6(k 2 - (u - u + (q)) 2 ) + 

+Z4q)(3_(u + k)(l - h(u.(q)))h(uj + u;_(q))9(k 2 - (u; + cu_(g)) 2 )] } , 

(28) 

where the trivial azimuthal integration has been performed and only the inte- 
gration over the polar angle 9 between k and p is left out. Finally the cut-cut 
contribution is given by: 

+1 oo +fc 

a cc (uj,p) = ^—j J dcosO J dkk 2 J duj\ n{uji)n{uj — uj 1 )9{q 2 — (uj — u)\) 2 ) 

-1 o -k 

(e^ - 1) {(1 + k • q) \j3 + {u u k)(3 + {u -u u q) + k)[3_(u - u u q)} + 

+(1 - k • q) {p+iuu h)P-(u - u u q) + P-(u u k)P + {u - u u q)}} . 

(29) 

For the sake of comparison we give below the expression of the free MSF at 
finite momentum, reported also in [11], which can be obtained, for massive 
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quarks, inserting into Eq. (19) the free quark spectral function 



p F (K) = (If+m) — [5(k - e k ) - 5(k + e k )\ 



(30) 



One obtains: 



a£Ju;,p) = N c N f J 



d 3 k 



■x 



\ ( k * Q TYl^ \ 

x { ( 1 + ± + ) (l-n(e k )-n(e q ))[S(u;-(e k +e q )) - 5(u + (e k +e q ))] + 



+ 1 



k • q m 2 \ 



n(e q )-n{e k ))[5{uj-{e k -e q )) - 5(u>+(e k -e q ))] \ , (31) 



where e k = y/k 2 + m 2 and K = (k , k). 

The delta functions can be exploited to perform the integration over the mo- 
mentum k. The general formula for the free spectral function for an arbitrary 
quark mass m and number of flavors Nf is then: 



<ee(^P) = ^(- 2 -p 2 )x 

x | 6{uj 2 - p 2 - 4m 2 ) 



/ i 4m 2 2 ^ 

CJ 2 — £> 2 



+ 6{p 2 -uo 2 ) — B \ , (32) 



where (5 — 1/T and 



A = log 1 + e 



-log 1 + e 



-!I«-Pa/1-J?3 



, (33) 



£? = log 



v y 



-logii + ZH-v^ 



(34) 



For m = we have: 



TV A/* 

r P s I,, c Si, .2 J ,2 ^ 



87T 2 



cu 2 -p 2 ) #(w 2 -p 2 ) 



1+ V 

p/3 . 



(35) 



In the limit where p — > Eq. (32) becomes: 



. N c N fn , 9 9 / 4m 2 / 2 

C(^0) = ^(^-4mVl/l - — (1 - -5- 



e~ + 1 



(36) 



in agreement with what quoted in [1,2]. 
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4-2 Finite momentum spectral function: numerical results 

In this Section we present our numerical results for the finite momentum HTL 
MSF in the pseudoscalar channel. We remind the reader that the HTL approx- 
imation has been shown to provide good results for different thermodynamical 
observables (namely in agreement with the ones given by lattice calculations) 
starting from T ps 2.5T c [26,27,28]. Nevertheless in the following, for the sake 
of completeness, we will often present results also referring to lower tempera- 
tures. 




10 20 

CO [fm 4 ] 



Fig. 5. The zero momentum HTL and free MSF (divided by uj 2 ) as a function of oj 
for two different temperatures. Also shown is the asymptotic high energy plateau. 

In order to have a reference benchmark, we start by considering the zero 
momentum case. In Fig. 5 we plot the spectral function divided by u 2 , in order 
to get rid of the uninteresting asymptotic high energy growth. Two different 
temperatures are considered. One can appreciate the dramatic difference in 
the behavior at low energy between the free curves (which vanish at uj = 0) 
and the HTL ones (which diverge for uj — > 0). One can also recognize the Van 
Hove (V.H.) singularities in the HTL curves arising from a divergence in the 
density of states due to the minimum in the plasmino dispersion relation. All 
the curves approach for large values of u the same high energy, temperature 
independent, plateau. 

We now investigate how things change at finite spatial momentum. In Figs. 6 
and 7 we plot the pseudoscalar MSF (divided by uj 2 ) for p = 1 fm -1 at different 
temperatures. The following general features concern all the considered cases: 
the (massless) non- interacting result vanishes on the light-cone (i.e. for uj — p), 
at variance with the HTL curves, which stays finite there. On the other hand 
when p is finite both the free and the interacting curves diverge as uj — > 0. 
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Fig. 6. The finite momentum HTL and free MSF (divided by uj 2 ) as a function of 
uj at p = 1 fm -1 for two different temperatures. Also shown is the asymptotic high 
energy plateau. 




CO [fm 4 ] 

Fig. 7. The finite momentum HTL and free MSF (divided by uj 2 ) as a function of 
uj at p = 1 fm -1 for two different temperatures. Also shown is the asymptotic high 
energy plateau. 

Furthermore in the HTL case the V.H. singularities are washed-out by the 
angular integration in Eq. (25). Finally both the free and the HTL finite 
momentum MSF approach the asymptotic plateau for large values of uj. 

In Figs. 8 and 9 we display the MSF (divided by uj 2 ) at T = 2T C for a range 
of momenta from p = to p = 10 fm -1 . We can thus follow its momentum 
evolution. 
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0.2 
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asymptotic 
p = HTL 
p = free 
p = 1 HTL 
p = 1 free 



oo [frn 1 ] 

Fig. 8. The HTL and free MSF (divided by uj 2 ) as a function of lo at p = and 
p = 1 fm _1 for T = 2T C . Also shown is the asymptotic high energy plateau. 



o.i 



3 



asymptotic 
p = HTL 
p = free 
p = 4 HTL 
p = 4 free 
p= 10 HTL 
p = 10 free 




co [fin 1 ] 

Fig. 9. The HTL and free MSF (divided by u 2 ) as a function of uj at T = 2T C for 
different values of the spatial momentum (in frn -1 ). Also shown is the asymptotic 
high energy plateau. 

In Figs. 10, 11 and 12 we display, for a range of temperatures from T = T c to 
T = 10T C , the finite momentum HTL MSF for p = 1,4, 10 fin" 1 respectively. 

Furthermore in Figs. 13, 14 and 15 we display how, at a given temperature, 
an increasing value of the spatial momentum affects the behavior of the MSF. 

In Figs. 16 and 17 we show the different contributions [pole-pole (pp), pole- 
cut (pc) and cut-cut (cc)] to the finite momentum MSF for two different 
temperatures. It clearly appears that, for large enough frequencies, the pp term 
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<B/T 

Fig. 10. The finite momentum HTL MSF (in units of T 2 ) as a function of uj/T for 
p = 1 fm -1 at different temperatures. 
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Fig. 11. The same as in Fig. 10 for p = 4 fm 
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Fig. 12. The same as in Fig. 10 for p = 10 fm 1 . 
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Fig. 13. The HTL MSF (in units of T 2 ) as a function of uj/T for different values of 
the spatial momentum. The plot refers to T = 2T C . 




co IT 



Fig. 14. The same as in Fig. 13 for T = 4T C . 




p = 10 free 



co/T 

Fig. 15. The same as in Fig. 13 for T = 10T C . Also shown are the (massless) 
non-interacting curves, which vanish on the light-cone. 
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Fig. 16. The different contributions (pp,pc and cc) to the finite momentum pseu- 
doscalar MSF for p = 4 fm _1 and p = 10 fm _1 , at T = 2T C . 
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Fig. 17. The same as in Fig. 16, at T = 4T C . 
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plays the dominant role: hence a closer inspection to the latter is mandatory. 




Fig. 18. The major processes contributing to the pole-pole term of the HTL MSF. 
The plot is given for the case p = 1 fm _ l and T = 2T C . 





Fig. 19. The other contributions to the pole-pole HTL MSF for p = 1 fm 1 and 
T = 2T C . 

Accordingly in Figs. 18 and 19 we display how the different processes quoted 
in Eq. (26) contribute, at a given value of the temperature and of the spatial 
momentum, to the pole-pole term of the HTL MSF. In agreement with charge 
conjugation symmetry, it turns out that the processes 2-3, 7-10, 8-12, 9-11 and 
5-6 give the same numerical contribution. 

In analogy to the zero- momentum case studied in [1,2], it turns out that the 
dominant process at low energy (at zero momentum in fact it's the only one) 
is the decay q + — > q- + M (together with its charge conjugate one); then the 
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pp term displays a wide gap for an intermediate range of energy till when 
the plasmino-antiplasmino annihilation starts contributing. Such a process 
initially gives a quite large contribution due to the large density of states; 
then it decreases rapidly with the energy because of the very small value of 
the plasmino residue. 

From Fig. 18 it appears that, for large enough frequencies, the dominant role is 
played by the quark-antiquark annihilation. Such a process starts contributing 
for uj larger than a threshold depending on the thermal gap mass m q acquired 
by the quarks in the thermal bath. 

On the other hand the processes 9-11 and 5-6 (q_ — > g_ + M and g_ — > q + + M 
together with their charge conjugate ones) turn out to be totally negligible, 
due to the low value of the plasmino residue and to the very small available 
phase space. 



4-3 Finite momentum temporal correlator 



In this Subsection we present our findings for the finite momentum temporal 
correlator defined in Eq. (5). 

In order to assess the impact of the interaction it is convenient to consider the 
ratio 

dw ct htl (w,p)K(w,t) 

. (37) 



GnTL(-iT : p) 



o 



Gfree(-ir,p) +f 

J du a iTCC (uj,p)K(uj,T) 



Such a ratio turns out to be finite for every value of r. In particular in the 
previous Subsection we showed how the non-interacting and the HTL MSF 
display the same high energy asymptotic behavior, both at zero and at finite 
spatial momentum. This implies that the above ratio approaches 1 for r — > 
(or ($). In fact, due to the structure of the thermal kernel given in Eq. (5), in 
this limit the correlator G{r) is dominated by the high energy behavior of the 
MSF. 

In Fig. 20 we start by considering the zero momentum case (already studied 
in [1]), showing the ratio Gwt~l{t) / G{ ree {r) for a range of temperatures from 
T = T c toT = 10T C . 

Then we move to the finite momentum case. In Figs. 21 and 22 the spatial 
momentum is kept fixed to the values p — 1 fm -1 and p = 4 fm -1 , respectively, 
and the temperature is let to change in the range from T = T c to T = 10T C . 
As the temperature increases the ratio GHTL( r )/G ? f rcc ('r) moves closer to 1, 
reflecting the running of the coupling, and the bump at r = (3/2 is smeared 
out. 
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x/p 



Fig. 20. The ratio Ghtl(t) / 'Gfree( T ) f° r different temperatures at p = fm 




x/(3 



Fig. 21. The ratio Gutl(t) / Gf iee (r) for different temperatures at p = 1 f m 1 . 

In Fig. 23 the temperature is kept fixed at the values T = 2T C , T = 4T C 
and T = 10T C , respectively, and we display how the ratio G , HTL('r)/G ! f r ee(T) 
is modified by changing the spatial momentum. The finite momentum effects 
are less and less important as the temperature increases, since the ratio p/T 
gets smaller and smaller. 

The role of the not vanishing spatial momentum in the interacting theory can 
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Fig. 22. The ratio Ghtl(t) / G{iee( T ) f° r different temperatures at p = 4 fm 1 . 
be also made visible by considering the ratio 

+oo 

/ duj a HTL (u,p)K(u,T) 
<?htl Ht,p) i 

G HTL (-ir,0) +~ ' W 

J du an TL (u,0)K(u,T) 



which is displayed in Figs. 24 and 25. It is clearly apparent the substantial 
quenching of the HTL meson propagator (with respect to the zero momentum 
result) occurring, for large spatial momenta, around r = (3/2. 



5 Conclusions 



We have examined the impact of finite values of the spatial momentum on the 
spectral density and on the temporal correlation function of a pseudoscalar 
meson for temperatures above the deconfinement phase transition and zero 
chemical potential. This amounts to study the properties of an excitation car- 
rying the quantum number of a meson and propagating (i.e. being not at rest) 
in the heat-bath frame. 

In our treatment we employed HTL resummed quark propagators. This has 
allowed us to perform many (but, of course, not all!) calculations analytically 
and to identify the different physical processes contributing to the MSFs and 
to the temporal correlators, a task impossible to achieve, of course, when the 
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Fig. 23. The ratio Ghtl(t) / Gf rcc (T) for different values of the spatial momentum 
at T = 2T C , T = 4T C and T = 10T C . 
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Fig. 24. The ratio Ghtl (t, p) /Ghtl (t, p = 0) for different temperatures at p = 1 
fm _1 (left) and p = 4 fm^ 1 (right). 




Fig. 25. The ratio Ghtl (t, p) /Ghtl (t, p = 0) for different values of the spatial 
momentum (measured in fm _1 ) for T = 2T C (left) and T = 4T C (right). 

same quantities are evaluated on the lattice (apart from the difficulties related 
to the MEM technique to extract a spectral density from an Euclidean corre- 
lator). 

In our investigation we have explored a range of spatial momenta from p = 
fm _1 to p — 10 fm _1 and of temperatures from T = T c (for the sake of 
completeness) to T = 10T C , although the HTL approximation is known to 
hold starting from temperatures of the order of T ~ 2.5T C . Clearly no intrin- 
sic difficulty exists in extending our calculations below such a temperature, 
but the increase of the gauge coupling g(T) as T moves downward makes the 
separation of scales (hard and soft) which lies at the basis of the HTL approx- 
imation no longer warranted. 

In summarizing our results we start by reminding the reader that, at zero 
momentum, the main features of the MSF are the presence of the Van Hove 
singularities (reflecting a divergence in the density of states) and the radi- 
cally different low energy behavior of the HTL predictions with respect to the 
free case. Such a contrast results particularly visible in the plot of a(u, 0)/c<j 2 , 
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where, for lo — > the HTL curve diverges while the free result vanishes. 
Moving to finite spatial momentum we first notice that many more terms, with 
respect to the zero momentum case, contribute to the MSF, reflecting the dif- 
ferent physical processes which may occur and which have been discussed in 
detail in the text. 

As a result of the calculation it turns out that the Van Hove singularities, 
so prominent in the zero momentum case, are smoothed out by the angular 
integration when p ^ 0. 

Another finding worth to be pointed out is that while the free MSFs vanish 
on the light-cone, the HTL curves stay finite. 

Finally the impact of a finite value of the spatial momentum has been investi- 
gated also for the HTL temporal correlator G(r,p), for which we plotted the 
ratios with respect to the free and the zero momentum results. 
In particular over the whole range of temperatures and momenta examined 
we found that the ratio of the HTL result with respect to the non-interacting 
correlator differs from one just for a few percent. 

Hence, although we do not present here the Fourier transform (along the z- 
axis) of the finite momentum correlator (integrated over all the values of r), 
which indeed represent a major numerical effort (work is in progress in this 
direction) we do not expect dramatic differences with respect to the free re- 
sult. 

We hope our work can provide a complementary and independent approach 
to the studies of MSFs and temporal correlators performed on the lattice. 
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A Adjustment of the parameters 

For the Nf = 2 transition temperature, following [29], we adopt the value 
T c = 202 MeV. The ratio T c / 'A m = 0.7721 for the N f = 2 case is taken from 
Refs. [30,31,32]. 

The running of the gauge coupling is given by the two-loop perturbative beta- 
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function, leading to the expression: 



g-\T) = 26 log 




+ — log \ 2 log 




} 



(A.l) 



where 



16vr 2 
1 



11 - 



4) 



(167T 2 ) 2 



( 



102 -38^ 
3 



) 



(A.2) 



The choice of the renormalization scale /i should reflect the typical momentum 
exchanged by the particles which, in an ultra-relativistic plasma, is of order T. 
For what concerns the precise numerical coefficient one should in principle let 
it vary within a reasonable range in order to get an estimate on the theoretical 
uncertainty of the calculation. Here, for the sake of simplicity, we simply adopt 
the choice /i = 1.1427rT suggested in [29]. 

We note that the running coupling g(T) defined in Eq. (A.l) enters our results 
only through the thermal gap mass of the quark m q = g(T)T/\/6. 
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